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along the treatment pattern to produce the desired overlapping effect. 



(57) Abstract: Scanning spot 
locations are generated for ablating 
tissue using a scanning laser beam 
(104) over a treatment region by 
fitting a target function representing 
a desired lens profile of scanning 
spots in a particular treatment pattern. 
Symmetry effects are utilized to 
simplify the process for determining 
the scanning spot locations. In some 
embodiments, the basis function is a 
two-dimensional function representing 
a two-dimensional section of a 
three-dimensional treatment profile 
which has symmetry a with respect to 
the two-dimensional section extending 
along the treatment pattern. For 
example, the treatment pattern is 
generally straight for myopic and 
hyperopic cylinders, and. is generally 
circular for myopia and hyperopia. The 
target function and the basis function 
may be discrete for implementation in" 
a software algorithm (20), and be fitted 
using a least square fit. The fit produces 
ablation depths for discrete scanning 
spots which are used to calculate' the 
number, of pulses at each reference 
position along the two-dimensional 
section. The pulses are distributed 



WO 01/67978 



PCTAJS01/08337 



GENERATING SCANNING SPOT LOCATIONS FOR LASER EYE 

SURGERY 

This application is based on and claims the benefit of U.S. Provisional 
Patent Application No. 60/189,633, filed March 14, 2000, the entire disclosure of which 
is incorporated herein by reference. 

BACKGROUND OF THE INVENTION 
The present invention relates generally to tissue removal techniques and, 
more particularly, to generating Ideations for scanning with a laser to achieve a desired 
ablation profile for correction of errors in vision during laser eye surgery. 

A scanning system has the ability to trace out an arbitrary pattern with a 

small low energy spot. In most cases, a small spot equates to finer scanning details but at 
the expense of requiring more pulses to remove a given volume. The notion that a small 
spot will give a better fit than a large spot is generally true for arbitrary spot and ablation 
shapes, but the spot shapes can also affect the fit. For instance, trying to fit round disks 
into a square shape will result in a residue. The emphasis on getting a good fit should be 
on choosing a good balance of spot geometry and size. 

Small spot scanners have their problems. A small spot will have smaller 
coverage per shot, thereby requiring more pulses to remove a given ablation volume 
(inversely proportion to the area of the spot size). A larger spot will require a 
substantially smaller number of pulses but at the expense of resolution. Understanding 
how treatment varies with spot overlap will ease our ability to create ablation patterns. 

Present ablation algorithms follow a removal or subtractive process. An 
ablation spot is placed on a location along the two-dimensional corneal surface. This spot 
ablates a volume of tissue from the surface to produce a crater on this surface. Another 
spot is applied at another location along the two-dimensional corneal surface and another 
crater is produced. The challenge lies in placing the craters in the correct locations such 
that adding up the overlapping craters will produce the desired surface without producing 
undesirable residues and requiring excessive processing time. 



WO 01/67978 PCT/USOl/08337 

SUMMARY OF THE INVENTION 
; The present invention relates .to tissue ablation utilized in, for instance, 

corneal sculpting, In general, embodiments of the invention determine a treatment table 
containing scanning spot locations and characteristics (e.g., size, shape, and depth) for 
5 directing overlapping scanning spots of a laser beam to achieve a target ablation profile. 
Some embodiments of the invention use symmetry effects (e.g., axisymmetry or bilateral 
symmetry) to simplify the overlap analysis to determine scanning spot locations for 
directing the laser beam to achieve a desired ablation profile. In specific embodiments, ' 
the spot center locations define a series of linear overlapping ablation paths, with the 

10 shape of the ablation paths preferably being selected in response to a desired change in an 
optical characteristic of the eye. The scanning spots can be uniform or variable in size, 
shape, and/or depth. Advantageously, these paths can produce ablation treatment profile 
elements which can be combined, adjusted, and positioned over the two-dimensional 
corneal surface to, in combination, produce the desire three-dimensional sculpting. The 

1 5 ablation profiles of the overlapping paths can be represented by basis functions, which 
may include an array of ablation profile data or can be expressed analytically for certain 
profiles. 

An aspect of the present invention is directed to a method of generating a 
! treatment table for ablating tissue using a scanning laser beam for generating scanning 

20 spots over a treatment region larger in area than the scanning spots. The method 

comprises providing a target function or lens function representing a desired lens profile 
for ablating the tissue by scanning spots of the laser beam on the tissue. A basis function 
represents a treatment profile produced by the overlapping scanning spots in a treatment 
pattern. The target function is fitted with the basis function to obtain the treatment table 

25 including scanning spot locations and characteristics for the overlapping scanning spots 
of the laser beam. In a specific embodiment, the scanning spot locations are randomized 
to produce a random scanning order. The scanning spot characteristics of a scanning spot 
may include scanning spot size, shape, and depth of the scanning spot at a specific 
scanning spot location. 

30 Any desirable treatment pattern for scanning with overlapping scanning 

spots of the laser beam may be specified. In some embodiments, the overlapping 
scanning spots define a treatment pattern in the form of, e.g., a linear path, and the basis . 
function is a two-dimensional function representing a two-dimensional section of a three- 
dimensional treatment profile which has symmetry with respect to the two-dimensional 

: ' 2 
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section extending along the treatment pattern. For example, the treatment pattern is 
generally straight for myopic and hyperopic cylinders, arid is generally circular for 
myopia and hyperopia. In a specific embodiment, the lens function represents an ablation 
depth as a function of a distance from an optical axis of a cornea. The basis functions of a 
series of offset patterns are combined, the profiles of laterally adjacent paths often 
overlapping to provide a smooth treated surface. 

In some embodiments, the spot size and shape are generally fixed. In 
other embodiments, the spot size and shape are variable. The ability to achieve the 
desired ablation profile using fixed spot size and/or shape may allow the use of a more 
simplistic laser source and simplify the ablation process. 

In specific embodiments, fitting the lens function with the basis function 
includes fitting at //discrete evaluation points. The basis function includes M discrete , 
basis functions representing M overlapping scanning spots. The M discrete basis 
functions may represent M overlapping scanning spots across a treatment zone length 
representing the length across a generally two-dimensional section which is oriented 
normal across a generally straight treatment pattern or which is oriented radially across a 
generally circular treatment pattern. 

For a treatment profile having a generally uniform two-dimensional 
section oriented normal across a generally straight treatment partem, the discrete basis 
functions represent the two-dimensional section as 

X t {xj) = y,(xj) = ^{sllf ^Xj -.x 0l .) 2 ; or, 

for a treatment profile having a generally uniform two-dimensional section oriented 
radially across a generally circular treatment pattern, the discrete basis functions represent 
the two-dimensional section as 



= 3(*,) = cos- 



V+x 0 2 ,-(s/2)^ 
2 ' x oi ' x j 



where 



s is the diameter of the scanning spot; 

j= 1.....N; 

Xj is a reference x-coordinate for the two-dimensional section measured 
from an optical axis of the cornea of a j* evaluation point for the center of the scanning 
spot; 

x 0 i is an x-coordinate for a center of an i th scanning spot; 
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(x Qi - sIT) < xj < (x 0i + j/2); 

yi(xj) is a depth of the i ,h basis function for the generally straight treatment 

pattern; and 

di(xj) is a coverage angle of the i th basis function for the generally circular 
5 treatment pattern. 

In a specific embodiment, fitting the lens function with the basis function 
comprises solving the following equation for coefficients a,- representing treatment depth 

for the i th scanning spot: 

•• m 

/(*,)=2a f X ( .(x,) 

i=i ' 

10 where 

f(xj) is the lens function; and 

. i= 1,...;M. ' 

Fitting the lens function and the basis function may include specifying a 
deviation for each of the N discrete evaluation points. The method may include refitting 

15 the lens function with the basis function by varying the deviations to iterate for an 

acceptable fit or a best fit. The method may further include refitting the lens function 
with the basis function by varying the number of scanning spots M to iterate for a best fit. 
The method may also include refitting the lens function with the basis function by varying 
the size of the scanning spot and/or the number of scanning pulses at a scanning spot 

20 location to iterate for a best fit. A merit function may be defined to represent an error of 
fit between the lens function and the basis function. The method includes minimizing the 
merit function to achieve the best fit. In some embodiments, the merit function and the 
total number of scanning spots in the treatment table are both minimized to achieve the 
best fit With the least number of scanning spots. 

25 Another aspect of the invention is directed to a method of generating 

scanning spot locations for ablating tissue using a scanning laser beam for generating 
scanning spots over a treatment region larger in area than the scanning spots. The method 
comprises providing a lens function representing a desired lens profile for ablating the 
tissue by scanning spots of the laser beam on the tissue. A basis function represents a 

30 treatment profile produced by the overlapping scanning spots along a treatment path. The 
basis function represents a section oriented across the treatment path. The lens function is 
fitted with the basis function to obtain the scanning spot locations for the overlapping 
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scanning spots of the laser beam. In specific embodiments, the treatment profile is 
symmetrical with respect to an axis of symmetry. 

In accordance with another aspect of the invention, a method for fitting a 
three-dimensional target profile comprises providing a two-dimensional basis function 
including overlapping portions to represent a three-dimensional profile which has 
symmetry with respect to a two-dimensional section extending along a treatment pattern. 
The three-dimensional target profile is fitted with the two-dimensional basis function to 
obtain a distribution of the overlapping portions. 

In accordance with yet another aspect of the present invention, a system 
for ablating tissue comprises a laser for generating a laser beam, and a delivery device for 
delivering the laser beam to a tissue. A controller is configured to control the laser and 
the delivery device. A memory is coupled to the controller, and comprises a computer- 
readable medium having a computer-readable. program embodied therein for directing 
operation of the system. The computer-readable program includes a first set of 
instructions for generating a treatment table for ablating the tissue over a treatment region 
larger in area than the spot size of the laser beam to achieve a desired lens profile for 
ablating the tissue, a second set of instructions for controlling the laser to generate the 
laser beam, and a third set of instructions for controlling the delivery device to deliver the 
laser beam to the tissue at the scanning spot locations. 

In specific embodiments, the first set of instructions includes a first subset 
of instructions for providing a target function representing the desired lens profile for 
ablating the tissue by scanning spots of the laser beam on the tissue, a second subset of 
instructions for providing a basis function representing a treatment profile produced by 
- the overlapping scanning spots in a treatment pattern, and a third subset of instructions for 
fitting the target function with the basis function to obtain the treatment table including 
scanning spot locations and characteristics for the overlapping scanning spots of the laser 
beam. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a top projection of a set of overlapping square pulses illustrating a 
straight pattern; 

Fig. 2 is a cross sectional projection of the overlapping square pulses of 
Fig. 1 along the central meridian illustrating the cross-sectional treatment profile; 

Fig. 3 shows a top projection of a treatment profile of a round spot trench; 
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Fig. 4a is a schematic view of a circle rotated about the z-axis; 

Fig. 4b shows a cross-sectional projection of an arbitrary shape; 

Fig. 4c is a top projection of a set of overlapping circular pulses illustrating 

a ring pattern; 

5 Fig. 5 shows the asymmetric treatment profiles of the overlapping circular 

pulses in the ring pattern of Fig. 6 taken along a radial section of the ring pattern; 

Fig, 6 is a one-dimensional, simplified representation of a target profile 
fitting equation; . ; 

Fig. 7 is a flow diagram illustrating a simulated annealing method for 
10 solving a target profile fitting equation; . 

Fig. 8A is a block diagram of the laser system illustrating an embodiment 
of the present invention; 

Fig. 8B is a block diagram of the control structure of a computer program 
according to a specific embodiment; . , 

15 . Fig. 9 shows the result of fitting a myopic lens equation using 27 equally 

spaced basis functions; 

Fig. 10 shows the result of an extended zone fit of a myopic lens equation 
with 47 equally spaced basis functions; ' ■ . 

Fig. 1 1 is an interferometer image of a -4 diopter myopic lens ablation 
20 profile using the simulation output from Fig. 10; 

Fig. 12 shows ablation profile changes due to spot size variations; ' 
Fig. 1 3 is an interferometer image of a hyperopic lens ablation profile; 
Figs. 14A and 14B show extended hyperopic lens fit with two different 

spot sizes; 

25 Fig. 15 is an interferometer image of a -4D hyperopic cylinder ablation 

profile; 

Figs. 16a- 16c are graphs illustrating the results Of the simulated annealing 
process at different stages of the iteration for fitting an elliptical profile; 

Figs. 17a and 17b show contour plots of the elliptical target profile and the 
30 calculated profile; and 

Fig. 18 shows a -4 diopter.ellipticai ablation profile based on the solution 
obtained from the simulated annealing process of Figs. 16a- 16c. " 
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DESCRIPTION OF THE SPECIFIC EMBODIMENTS 

I. : Generation of Scanning Locations and Depths 

A. Treatment Space 

Treatment space or profile is the treated form (or profile) created by a 

5 given spatial pulse sequence taking into account its overlap characteristics. 

1. Rectangular Spot Shap e in Straight Pattern 

As an introductory example, a simple square spot with a treatment 

intensity of approximately .25 p.m per pulse is used. The overlap is set to 33%. A 1 mm 

square spot centered at position (in mm) [0, 0] will have the next spot at [0, 0.33]. The 

10 area the two pulses have in common is 0.67 mm 2 . The offset (step size) from the first 

pulse is 0.33 mm. Continuing this sequence for a total of N pulses will produce a treated 

path of width 1 mm, and length 0.33*N mm. The overlap is 33% so the last pulse to 

coincide with the first pulse is the third pulse. The fourth pulse's edge is immediately 

adjacent to the first pulse. Assuming that the positioning is perfect, they will share a 

15 common edge. Looking in from a meridian cross section, the slope is (depth per 

pulse)/0.33 mm. All positions after the third pulse will have a depth of (3* 'depth per 

pulse). The slope also occurs at the trailing end. The cut profile from the ends has the 

box shape, with a depth of (3 *depth per pulse) and a width of 1 mm. The slope is: 

Slope = ^Pih per pulse 
step size 

20 The treatment profile of Figs, l and 2 is generated by scanning 

overlapping square spots in a straight treatment pattern. As shown in Fig. 1 , the shape of 
the trench is rectangular in shape. The stacking of pulses has the first interval with a 
single pulse depth, thfc second interval with a double pulse depth, and the third interval 
with a triple pulse depth, as seen in Fig. 2. The depth of the trench remains the same after 

25 the third pulse because the overlap is 1/3 of the spot, i.e., the fourth pulse does not 
overlap the first. 

The treatment profile generated by overlapping rectangular spots is 
relatively simplistic and easy to represent mathematically due to the rectangular spot 
shape as well as the straight-line pattern. The complexity of the mathematical 
30 representation increases for round spots, particularly if the spots overlap in a nonlinear 
pattern. 

2. Round Spot Shape in Straight Pattern 
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Fig^ 3 shows a round spot having a diameter or spot size of s. x is a spatial 
coordinate from the origin to a point on the trench produced by the round spot (which 
may be discrete). x 0 is the position of a basis function or the center of the ablation trench 
produced by the round spot, and 0 is the coverage angle. A basis function is a 
5 mathematical representation of the treatment profile produced, for instance, by 

overlapping scanning spots in a specified treatment pattern. Here, each round spot is 
assumed to produee an ablation with a uniform ablation depth. Adjustments may be 
made for nonuniform energy distribution profiles, "central islands," and other local 
ablation depth variations, and the like. 

10 Overlapping a round spot produces similar characteristics as overlapping a 

.^j rectangular spot. Repeating the samefexercise for overlapping the round spot in a straight 
pattern as done with the square spot, the side profile will have the same look as the square 
spot. The leading and trailing slopes will Have the same slopes as those of the square 
spot. The profiles along the ends are different. It is semicircular in shape, or more 

15 precisely, semi-elliptical in shape, with a maximum depth of 

* ^ spot size , ■ : 

Depths =— depth per pulse. 

step size 

The slope is the same as described in the square spot. The round spot has 
a different edge pattern. Without a better.descriptive phrase to describe it, it has a 
"cloud" effect to its edge. 
20 The basis function for a round spot In a straight pattern is: 

Yi 00 - VC^POt size/2) 2 - (x - x oi f 

with 

(x 0 i - spot size /2)<x< (x 0 i + spot size / 2); 

where 

25 / is the running index of the basis functions; 

xoi is the position of the i th basis function (i.e., the center of the i th spot) 
measured from the origin; 

yi(x) is the trench depth of the i th basis function; and 

x is the spatial coordinate from the origin (e.g., optical axis of a cornea) to 

30 a point on the trench. 
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Fitting the function in a computer algorithm requires the coordinate system 
to be discrete, Definingy as the running index for the x spatial dimension, the trench 
depth is the same as above but with x replaced by xj, 

yi( x j)= >J(jpotsize/2) 2 -(x. -x oi ) 2 

where j is the running index for the sampling or evaluation points. 
3. Round Spot Shape in Ring Pattern 

Generation of a ring treatment pattern makes use of axisymmetry. To 
illustrate how symmetry affects the overlap, consider a circular disk of radius r with its 
center b ' away from the origin as depicted in Fig. 4a. Recall from differential geometry 
that its normal vector (right-hand rule) gives the "direction of the surface". The disk's 
direction vector at anytime is in the x-y plane thus implying that the surface of the disk 
itself lies in the y-z plane. Revolving this disk about the origin in the direction of its 
normal vector (z-axis) will generate a bagel shaped volume. The bagel will have a radius 
of b ' and a thickness of 2r. This assumes that b '>r. 

Consider again, the same disk but only the top half (z>0). The bottom is 
cut off at the x-y plane. Rotating this half disk about the surface along its normal vector 
will generate a semi "flat bottom sliced bagel" ring with a diameter of 2b ' and thickness 
ofr. Integrating the volume as the half disk moves about the origin will produce the 
volume of the flat bottom bagel. This type of "differential volume" will allow us to 
produce a "bowl" shape (flat on one side, curved on the other) similar to the bottom half 
of a sliced bagel. 

One method for calculating volume is a 3-D version of Simpson's rule. 
Fig. 4b shows a cross section of a set of circular rings at a, 'x lt . . .,.x (n .j) with a constant 
thickness of Ax: The top curve is the desired target shape representing, for example, a 
lens volume. The columns are cross sections of "thick" rings, spaced xj apart. Summing 
a particular set of rings will produce the desired volume. Note that this is mathematically 
easy to do because the basis functions are orthogonal. The rings following Simpson's 
rule will need to be very thin to get a good approximation to the curve (i.e. the error is on 

the order of &(x)). 

Generalizing these approaches, symmetric objects can be generated in a 
"space" minus one dimension. For example, a symmetric 3-D object can be generated in. 
2-D space if symmetry is utilized. Generating these basis "rings" or treatment volumes 
from differently sized circular spots is the basis for one aspect of the present invention. 
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The thin lens volume is generated using flat bottom bagel shaped rings. 
Finding the different combinations of ring radii (b *) and thickness (r) to. generate the lens 
is done via, for example, a linear least square fitting routine. 

Scanning a continuous circular spot in a large ring pattern will produce a 
5 result that resembles a bagel with a half-circular cross section (like a sliced bagel). As the 
ring's position gets closer to the origin, the ring's cross section becomes very noncircular 
because the ablation spot overlap characteristics change as the ring's radius changes. 

The scan pattern for a round spot in a ring pattern will look somewhat like 
• a torus but with a radial dependence resulting in a bulge toward the origin. To understand 
10 this radial dependence, consider a spot moving along a ring as shown in Fig. 4c. The 
outer portion of the spot will have an angularly smaller coverage , than that of the inner 
, portion. That is, the overlap for the outer portion of the spot is less than the inner portion. 
The functional expression for the arc length for this spot as a function of the radial 
distance is an arccosine: 



or in terms of spot size or diameter, 

0. (x) = cos" 

where 



r x 1 + jcp. - (spot size/ 2) 2 ^ 

2 *o, * 



i is the running index of the basis functions; 
20 x is the spatial coordinate from the origin (e.g., optical axis of a cornea) to 

a point on the trench; and 

xoi is the position of the i, basis function (i.e., the center of the i th spot) 
measured from the origin; and 

0i(x) is the coverage angle of the i th basis function. 
25 Fig. 3 shows that the i ,h basis function is defined from x 0 j - s/2 <x <x 0 , + s/2 

and zero otherwise, x is continuous but actual implementation in an algorithm to solve 
for the scanning coordinates or locations will have discrete x. 

A discrete coordinate system facilitates fitting the function in a computer 
algorithm. Defining^' as the running index for the x spatial dimension, the angle of 
30 coverage is the same as above but with x replaced by xj, 
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' x. 2 - + xl -(spot size/ 2) 2 
*,(*,) = co.-M — 



where; is the running index for the sampling or evaluation points. 

A feature of these treatment "rings" is the asymmetric treatment space 
profile with respect to the radial direction from the origin as shown in Fig. 5. The profiles 
5 get more asymmetric as the pulse rings approach the origin, as illustrated by the increased 
asymmetry of profile A over profile B. This feature allows the fitting of a centrally 
maximum curved surface. The coefficient of each basis function will give the treatment 
depth. 

4. Treatment Matrix 
10 The above examples employ specific treatment profiles that can be 

expressed in analytical forms for achieving relatively simple types of target ablation 
profiles. Other treatment profiles may be used. A general treatment matrix includes data 
of treatment profiles produced by overlapping scanning spots of any desired shapes and 
sizes. . 

1 5 At one time it was thought that the laser-tissue interaction resulted in a flat 

ablation profile, so that each pulse could be" thought of as removing a disc of material, 
with a diameter equal to the laser iris width and a uniform depth. Empirical evidence, 
however, has shown that the laser-tissue interaction may result in a rather complicated 
ablation profile with a complex shape and characteristics that can vary significantly with 

20 the laser iris diameter. In general, the ablation profiles for different laser spot sizes can be 
digitized and arranged as a treatment matrix of ablation data to be used for numerically 
determining a treatment plan or table (i.e., ablation amplitudes for a treatment pattern of 
scanning spots) to achieve a target ablation profile. There will be a different set of data 
for each treatment medium including, for example, polymethylmethacrylate (PMMA), 

25 calibration plastic, human comea, porcine cornea, or the like. 

B. Target Profile Fitting 

Fitting the basis function to the target function representing the target 
ablation profile produces the scanning spot locations and ablation depths. In specific 
embodiments, the target function is represented by the lens equation as described below. 
30 The target profile fitting involves solving for a set of coefficients for the 

following target profile fitting matrix equation: 
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" /« = f> ; *,(*) 
. F = AX 

with the solution as an inversion 

a = x'f 

5 where , ' 

Xtfx) is the i th basis function; 

f(x) is a particular target profile, which may be expressed as a lens 
equation or represented by a target matrix of data; and 

a,- is the amplitude of the i* basis function. 

,10 F is the lens equation or target matrix, X are the basis functions, and A is the treatment 
plan or table containing the amplitudes for a particular treatment pattern (i.e., the 
scanning spot locations and characteristics such as size, shape, and depth). The fitting 
program calculates A such that the error (x 2, s) between F and X is a minimum. 

Some of the examples discussed above employ analytical basis functions. 

15 For instance, the basis functions.X are y/xj) for the straight pattern and (^fx^ for the ring 
pattern. The straight pattern is applicable in myopic and hyperopic cylinders, while the 
ring pattern is applicable in myopia and hyperopia. For photofherapeutic keratectomy 
(PTK), the treatment profile is essentially a centrally flat hyperopia ablation since PTK is 
a flat ablation- Thus, the basis function &,(xj) can be used for PTK (while the lens 

20 equation for PTK is a constant value). Of course, the basis function may be represented 
by a general treatment matrix instead of the analytical functions. Furthermore, the target 
profile/^ may be expressed as a lens equation which has been developed for certain 
cases such as myopia and hyperopia. Other target ablation profiles (e.g., an elliptical 
profile, a toric profile, a wavefront profile, or an arbitrary profile) may be expressed more 

25 generally as a target: matrix. 

Fig. 6 schematically illustrates in one dimensional simplified form the 
equation XA = F, where X is a basis function expressed in a general treatment matrix 
representing ablation profiles that may have relatively complex shapes, F is the target 
profile schematically shown as a lens, and A is the treatment plan or table to be solved. 

30 In this example, the treatment table is a pulse insfmction vector (PIV) representing the 
laser instruction count for each basis function. In the specific example shown, the 
treatment matrix represents ablation profiles for five laser pulses (i.e., five particular laser 
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pulse diameters or shapes at particular locations). The elements of the treatment table or 
PIV correspond to the numbers of pulses that can be combined to build an entire ablation 
having a shape that matches the target profile as close as possible. In that case, the 
elements of the treatment table or PIV are positive integers. The treatment table or PIV is 
used to operate the laser to achieve ablation that approximates the target profile. Note 
that a two-dimensional version of the equation (instead of the one-dimensional illustration 
shown in Fig. 6) will produce a three-dimensional model of the ablation profile. 

1 . Least Square Fitting Approach Using Analytical Len s Equations 
The difficulty in performing the inversion of the target profile fitting 
equation to obtain A is that X" 1 is not orthogonal . One useful method for solving this 
equation is to use a least square fitting routine. An example of a suitable method for 
fitting the basis function is found in Numerical Recipes in C, 2 nd Ed., W. H. Press, S.A. 
Teukolsky, W.T. Vetterling, B P. Flannery, Cambridge Univ. Press, pp. 671-675 (1992). 
The general model representation is 

/(x) = f>,*,(*) 

M 

where X\, . . ., X M are arbitrary functions of X, called basis functions, and M is the number 
of basis functions. The A^s need not be linear. The linearity is in respect to the linear 
dependence on the parameters a t . The basis functions are the v/'s for the straight 
treatment pattern or the 6? s for the ring treatment pattern. 
The merit function (x 2 ) is defined as 

r m l 2 

/(*;)-!>,-*,•(*,) 

2 _ M 



f(x) represents the target profile to be fitted with the basis functions. In' specific 
embodiments,/^ is the lens equation evaluated at*,, that is, f(xj). The a,-'s are the 
uncertainties or standard deviations of the Xfs. These are typically set to one but can be 
altered if a more precise fit for certain parts of the curve is desired (e.g., the central half 
can be set to 0.25 and the outer half can be set to 1). For instance, decreasing the value of 
a, for the central area and increasing it for the outlying area will force the fit to be more 
precise in the central region. The reduced % 2 is 

V 
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Where u is the number of degrees of freedom. 

To minimize jf, the derivative of with respect to parameters a* for all 
Afs is taken and set to zero. ". 

M 



1=1 o"i L '=i 

5 for i=l, M These are the normal equations of the least square problem. The two 

commonly used methods of obtaining the vector a is by LU decomposition and by Gauss- 
Jordon elimination. The Gauss-Jordon method will provide the variances of a,- as part of 
the solution. 

%l gives the relative "goodness" of fit. The better the fit, the closer 
1° zl approaches one (%l -» +1) . If zl is Iess than one, it means the original estimates of 
of 's are too large. 

The input parameters used by the fitting routine are: 

(1) The data generated by the lens equation, /fx^. . 

(2) The standard deviation pi of each data point. Since the data is 
15 generated from the lens equation, the standard deviation is arbitrarily set to one. This 

arbitrariness needs to be addressed because the standard deviation affects %* '. Consistency 
is desired because once the values for the standard deviation is determined, it should not 
be changed if x 2 is to be used to determine the relative goodness of fit. A wide variety of 
fit evaluation methods could alternatively be used within the scope of the invention 
20 (3) The basis functions in analytical form. 

2. Simulated Annealing Approach 

Another computational technique that may be used for solving the target 
profile fitting equation is simulated annealing. Such a technique is described, for 
instance, in Numerical Recipes in C, 2 nd Ed., W. H. Press, S.A. Teukolsky, W.T. 

25 Vetterling, B.P. Flannery, Cambridge Univ. Press, pp. 444-451 (1992). Simulated 
annealing is an evolutionary optimization method using gradual, randomly induced 
changes to find the best solution to mathematically intractable problems. 

Simulated annealing is a method used for minimizing (or maximizing) the 
parameters of a function. It is particularly suited to problems with very large, poorly 

30 behaved function spaces. Simulated annealing can be applied in the same way regardless 
of how many dimensions are present in the search space. It can be used to optimize any 
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conditions that can be expressed numerically, and it does not require a derivative. It will 
not be fooled by local minima in the search space. 

Simulated annealing is inspired by the observation that if a thermodynamic 
system is cooled slowly, it will end up very near its minimum energy state. For example, 
a piece of hot metal contains n molecules in rapidly changing orientations. If the metal is 
quenched in cold water the thermal energy rapidly leaves the metal and the molecules 
remain — "frozen forever"— in their last orientation. However, if the metal is allowed to 
cool gradually the molecules will continue changing orientation, but gradually favoring 
lower and lower energy states. Eventually all the molecules will be "frozen" into their 
lowest energy state. (In metals, the molecules align, forming crystals.) This is called 
annealing. 

Instead of a piece of metal, the treatment table or pulse instruction vector 
(PIV) will be annealed. The PIV may be expressed as a 1 by N array of integers, each of 
which describes the number of pulses at a particular location and given pulse size, shape, 
and/or depth. The PIV may be expressed as follows: 

tt of 1 mm pulses at 
position 1 

It of 1 mm pulses at 
position 2 

# of 1 mm pulses at 
position 3 



It of 6 mm pulses at 
position M-1 

tt of 6 mm pulses at 
i position M 

This vector can be very large. Some examples have upwards of 3000 elements. 

To "simulate" annealing, two physical concepts are applied to solving the 
target profile fitting equation. The first is the concept of energy state to be minimized, 
which is analogous to the error or merit function (x 2 ) defined as 

where f(x) represents the target profile to be fitted with the basis functions X{. In specific 
embodiments, /fx; is the target profile evaluated at x J} that is, f(xj). The numerical value 
(X 2 ) to be minimized is the sum of the squares of the difference at each point between the 
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target shape and the "test shape" or "working shape" during the iteration. This example 
employs a least square fit for minimizing the merit function. > 

The merit function may be modified to find a fit with more desirable 
optical qualities. For instance, a weighted merit function may be 



M . 



1=1 



f(x j )-a i X i 



-12 



The o,'s are the uncertainties or standard deviations of the X t 's. These are typically set to 
one but can be altered if a more precise fit for certain parts of the curve is desired (e.g., 
the central half can be set to 0.25 and the outer half can be set to 1). For instance, 
decreasing the value of cr, for the central area and increasing it for the outlying area will 

10 force the fit to be more precise in the central region. 

The second physical concept that is applied in the fitting approach is that 
of temperature, which is analogous to a measure of the random change applied to the 
treatment table or PIV (i.e., change in the number of laser pulses) before each successive 
evaluation. In this case, each "degree" of "temperature" change will cause one random 

15 fluctuation in the PIV. For each fluctuation one element of the PIV will be randomly 

selected. That element will then be either incremented or decremented, again based on a 
random number. In the case that the element is already zero and a decrement is chosen, 
nothing will happen to that element, since only positive integers are acceptable. 
Evaluation of the "energy state" (i.e., merit function) takes place during the "temperature" 

20 change iteration for the PIV until a sufficiently low value for the merit function is 
achieved. 

An example of implementing the simulated annealing procedure or 
algorithm is illustrated in Fig. 7. After the master application calls the algorithm (step 20) 
and the simulated annealing subroutine is initiated (step 22), the subroutine randomly 

25 modifies the treatment table or pulse interaction vector (PIV) according to the current 
temperature (i.e., size of the incremental change to an element of the PIV) to create the 
next vector PIV ' (step 25). An evaluation takes place in step 26 to determine whether 
PIV is better than PIV in reducing the merit function or "energy state." If PIV is better 
than PIV, then the vector is updated by setting it to PIV (step 28), and the iteration 

3 0 continues by following the loop back to step 24. If PI V is not better than PIV, then the 
subroutine determines whether the vector has remained constant for a preset number (n) 
of evaluations (step 30). If not, the PIV is not updated and the iteration continues by 
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returning to step 24. If n number of evaluations have taken place with no improvement of 
the PIV, then the subroutine determines in step 32 whether the temperature .(i.e., the size 
of the incremental change) has reached a preset minimum value. If not, the temperature 
is lowered by one degree in step 34 (i.e., the size of the incremental change to an element 
of the PIV is lowered by a preset amount) and the iteration continues by returning to step 
22. If the minimum temperature has been reached, then the annealing is complete (36) 
and the simulated annealing procedure is terminated, Therefore, the annealing process 
constitutes repeated evaluation of the PIV, at gradually decreasing temperatures. For a 
particular temperature the iteration cycle will repeat until a certain n number of 
evaluations have failed to show improvement. Once the solution meets some pre- 
determined "goodness of fit criteria" as represented by the merit function, the final 

solution is returned. 

For problems with radial symmetry, such as spheres or round PTK's, the 
problem is solved in only one dimension. For problems with bilateral symmetry, such as 
ellipses or cylinders, the problem becomes two dimensional. The method only needs to 
be used to solve one quarter of the surface and then the solution is applied oyer ail four 
quadrants. This saves considerable computation. For problems with no symmetry, such 
as a wavefront map, the method solves for the entire surface. 

Solving the target profile fitting equation numerically for a general two- 
dimensional. surface, particularly without the benefit of symmetry, can be very difficult 
and can require a great deal of computational time. The simulated annealing approach 
can generally produce the best solution in a reasonable amount of time, preferably in real 
time. 

The input parameters used by the simulated annealing routine are: 

(1) The target profile data, which can represent any arbitrary shape. 
The target shape can be arbitrary, allowing for unlimited ablation possibilities, including 
the map from a wavefront type device or any other topography system. As long as the 
target profile can be described mathematically or with an array of data, the algorithm can 
use it. Radial or bilateral symmetry is not required, although the presence of symmetry 
can make the computation faster. 

(2) The treatment matrix (i.e., basis function) which can represent any 
ablation sizes and shapes. The laser ahlation shape does not have to be flat or smooth, or 
have specific. shapes. As long as it can be measured and described mathematically or 
with an array of data" the simulated annealing process can be Used. A unique description 
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s ... 
is provided for every unique pulse shape or size to be used. A fixed spot-size laser would 

have, only one description, while a variable spot-laser could have as many as desired: 

(3) The standard deviation o\ of each data point if a weighted merit 
function is used. It is noted that any merit function can be Used. This allows the user to 
5 decide what aspect of the solution is most important. For example, extra "weighting" 
may be provided in the central portion of the lens, so mat the algorithm will recognize 
that the center portion of the lens is the most important part and will sacrifice accuracy 
around the edges to get the best possible fit in the center portion. Moreover, the merit 
function is not limited to just the shape of the surface. For instance^ the merit function 
10 can be linked to the number of pulses, so that the algorithm will favor solutions with 
fewer pulses. This helps to minimize the treatment time. In one example, the merit 
function is expressed as follows: 



• M 



-|2 



• X number _of ^pulses) + ^ offset) 



In this example, the merit function is just the surface fit function 
15 . multiplied by the number of pulses squared plus a constant offset. The offset lets the user 
tune the relative importance of the surface fit versus the pulse count. If the offset is large 
compared to the number of pulses, then the merit function will be only mildly influenced 
. by pulse count. The solution will tend to reduce the number of pulses, but only slightly. 
If the offset is comparable or small compared to the number of pulses then the merit. 
20 function will be significantly influenced. The solution will compromise some goodness 
of fit in order to use fewer pulses. It is possible to have solutions with differing pulses 
because a variety of pulse sizes are used (e.g., a single five millimeter pulses can remove 
as much tissue as twenty five pulses of one millimeter diameter). 

Other modification of the merit functions would include maximizing 
25 smoothness, (i.e. minimizing the first derivative of the surface error), meeting hardware 
constraints of the Laser delivery system (e.g., the Laser can only move so far in a given 
time), meeting biological constraints on the eye (e.g., pulses,rriight be distributed over 
space and time to avoid excess thermal build up.) 
C. Software Implementation Steps 
30 1. Least Square Fitting Approach Using Analytical Lens Equations 
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The following steps are provided to illustrate an example of software 
implementation of the fitting approach of generating scanning locations and ablation 
depths using analytical lens equations. 

Step 1. Evaluate Lens Equation. The lens equation is evaluated for N 

5 points (typically equally spaced) in the region 0< x <L where L is the treatment zone 

radius for a ring pattern and is the treatment zone length for a straight pattern. N is set to 

about 300 in specific embodiments. If it is desired to emphasize the central portion of the 

lens equation (e.g., for myopia and myopic cylinder), a "shift" is subtracted from the 

treatment zone length. 
10 For myopia and myopic cylinder, the lens equation is: 



where 



-x] +C 



0<xj <_(L - shift); and 
j = 0,1,..,,N-1. 

15 The constant C. is defined as 



where 

xj is the reference x-coordinate measured from the optical axis of a cornea 
of thej th evaluation point for the center of the ablation or scanning spot; 
20 s is the diameter of the ablation (i.e., spot size); 

Ri is the anterior radius, of curvature of the cornea in meters; 

R 2 is the final anterior radius of curvature of the cornea in meters; 

n = 1,377 is the index of refraction of the cornea; 

D is the lens power of the ablation in diopters; 
25 L is the treatment zone length representing the length across a generally 

uniform section which is oriented normal across a generally straight treatment pattern or 
is oriented radially across a generally circular treatment pattern; arid 

shift is the amount of emphasis shift {shift typically is about 0-0.2, and can 
be set to zero). ' 
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Note that the above lens equation applies for the fiill range of angles in 
myopia to provide an axi symmetric profile, but applies only at a particular angle in 
myopic cylinder to provide a constant profile along a straight pattern. 

For hyperopia and hyperopic -cylinder, the lens equation is: 



■ 71-1 + /?,/) V ' 



'{n-l + R^D J 
where 

0<xj<(L- shift); and 
j = 0,J,...,N-J. 

Note also mat this lens equation applies for all angles in hyperopia and at a particular 

10 angle in hyperopic cylinder. A.more detailed description of the lens:equatiori is provided 
in Charles R. Munnerlyn et al., "Photorefractive Keratectomy: A Technique for Laser 
Refractive Surgery," J. Cataract Refract. Surg. (Jan. 1988), which is incorporated herein 
- by reference in its entirety. 

For PTK, the lens equation is: 

15 f(xj)=d ' 

where d is. a constant depth. 

Step 2. Define basis rings at different radial positions for a ring 
pattern or define basis strips at different lateral positions for a straight pattern. A 
set of radial positions for basis rings or a set of lateral positions for basis strips are 

20 defined at *o/- In specific embodiments, the number of rings/Strips, M, may vary from, 
about 17 to about 57. The positions of the ring/strip's radial center are equally spaced 
points in the region 0< x < (L - spot size/2 + extended zone). The extended zone is for the 
taper or "feathering" of the outer circumference or edge; The extended zone typically 
ranges from about 0. 1 -0.5 mm, arid may be fixed at about 0.2 mm. The extended zone 

25 also helps in the fitting routine. M is the number of basis functions. The number should 
be large enough to give a smooth fit and not too large to create quantized errors. 
Typically, M is about 7-97. In specific embodiments, M starts at about 17 for -1 .5D or 
less, becomes about 27 between -1.5 and -2.5 D, and reaches about 47 for greater than - 
2.5 D. The i th basis function is 

30 xqi = [(treatment zone length ■ — spot size + extended zone) / M] * / 

wherein / = 1,2, M. 
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Step 3. Define the positions for lens equation evaluations. N is the 

number of evaluation points. The i th position for lens equation evaluation is 

xj = [treatment zone length — shift) / N] * j 

where7 = 0, 1, .... N-l. 

5 step 4. Assign a deviation, o> for each xj. The deviations will determine 

how tight the fit will be for each point. In a specific embodiment, the deviations are set to 
about 0.25 for the first 200 points as 0.25 and about 1.0 for the final 100 points. The 
deviation can be changed and optimized for different profiles and patterns. 

Step 5. Apply least square fit. The function to fit the basis functions to 
10 the lens equation (i.e., fit X(xj ) to flxj)) is 



where, for a ring pattern (myopia, hyperopia, and PTK), 

x / + x oi ~ ( s P ot size 1 2 ) 2 



X i (x J .) = 0 i (x J ) = cos- 1 



2-x 0i -Xj 



and, for a straight pattern (myopic cylinder and hyperopic cylinder), 
15 X,(x'j) = y t (xj) = ^(jpotsize/2) 2 -( Xj -x 0i ) 2 

Note that cos" 1 is valid only in the range (x 0 , - spot size / 2) < xj < (x 0i + spot size / 2) for 
the ring pattern. The index j is the running index for evaluation and fitting points. The 
index i is the running index for the number of basis rings or basis strips. 

The least square fit routine is an object that receives xj coordinates, the lens 
20 equation evaluated at these coordinates (f(xj) values), the basis functions Xfxj), and 

returns a row (or column) matrix A representing the coefficients, and a Once again, a 
variety of fit or optimization routines might alternatively be employed. 

Step 6. Evaluate closeness of fit. is recalculated if necessary, that is, 
if any changes are made to the returned matrix A. £ will give the goodness of fit. The 
25 difference between the lens equation and sum of the basis functions (standard deviation, 
etc.) can be evaluated to check for cusps. Cusps are small deviations with sharp 
transitions typically occurring at an integer multiple of the spot size. The cusps tend to be 
"stiff equations and are difficult to eliminate later. Cusps tend to be localized, occurring 
within 20 or so coordinate points, A moving average and standard deviation "window" of 
30 about 20 points can be done to test for cusps. Note that the effects near the edge of the 
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ablation will typically have a large standard deviation! If the deviation due to cusps is too 
large, different fitting parameters can be used to re-evaluate the closeness of the fit. The 
results should be checked for negative coefficients and their magnitude. Another fit with 
a different set of rings (usually a smaller number but may be a larger number) will need to 
5 be done if zeroing out the negative coefficients causes a large For small values, it 
suffices to set the negative coefficients to zero and re-evaluate £ and check the lens 
equation versus the solution. 

Step 7. Iterate for best solution. If there are 2 cusps (criteria to be 
determined) or if is too large (>100 for the initial settings of 0.25 and 1.0), or if the 
10 solution is not close to the lens equation (to be determined), refitting is performed using 
different fit parameters: 

(A) Begin with the number of basis functions, searching from 17-57 at 
the initial spot size (2 mm); 

(B) Next search with a different spot size (L75 mm, 1.5 mm) and 

15 repeat (A); 

(C) Next search with extended zone and treatment size such that the 
sum of extended zone and treatment dimensions is a constant (i.e., 6.2 mm) in 0.05 mm 
increments, and repeat (A) and (B); and 

(D) Next search with different o/s. This can get rather complicated. 
20 In one calculation, the initial 200 qj's were 0.25, and the final 100 qj's were 1.00. The 

initial 200 oj's were changed to 0.5 arid (A), (B), and (C) were repeated. Then the 200 
<^'s were changed to 199 and (A), (B), and (C) were repeated, and so on. 

Step 8. Determine actual pulse locations. The coefficients a t returned 
from me fit (matrix A) are used to calculate the number of pulses at each 
25 for a ring pattern or each lateral position for a strip pattern (i.e., basis functions), 

n* = a* / (depth per pulse) 
With the present fitting routine, if the coefficient is less than zero, it is typically set to 
zero. For myopia, the coefficients are rarely zero. 

For a ring pattern, the pulses are spaced evenly along the circumference, 
30 but with a random angle added to the first pulse. Dividing this number into 2ii will give 
the angular spacing between each pulse^ The first pulse of each ring is desirably ■• 
randomized to avoid a pattern effect.. " 
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Step 9. Randomize table. All pulses are grouped into one and 
randomized before treatment to avoid sequentially generating spatially adjacent treatment 
pulses. 

Geometrical considerations can be addressed during the fitting process for 
a ring pattern (as in myopia or hyperopia) or after the fitting process for a straight pattern 
(as in myopic and hyperopic cylinders). For reference, the equation for the treatment 
space profile for the ring pattern in the fitting step is repeated here: 

x lt + x ] ~ (spotsize 1 2) 2 



a i 0 i (x J .) = a i cos- i \ 



2 ■ x oi - Xj 



The at's are the basis function amplitudes returned by the fitting function 
10 (step 7). Each amplitude represents the height of a "ring". The length of the ring is 2tdc 0 u 
■ giving the total number of pulses as for the ring as follow: 

total pulses = : — : — 

spot size depth per pulse 

Because the basis functions have radial dependencies, it is easier to 
implement the normalization factor during the fit by using the normalized basis functions 
15 in the fitting step (step 5) as follows: 

f x li + x ) - ( s P ot size / 2 ) 2 > 
spot size • Xj j 

The a,'s returned from these basis functions are normalized. The total 
number of pulses per ring is , 

a , 

total pulses = 



a^.(x.) = ^^.(x,) = ^cos- 1 



depth per pulse 

20 The pulses are spread evenly along the 2n arc. According to a specific 

embodiment of randomizing the pulses for treatment in step 9, random phases are added 
to the starting points (the first pulse per radii). The random phases will reduce the 
periodicity caused by the "finiteness" of the pulse number. 

2. Simulated Annealing Approach 

25 The steps for implementing the method of generating scanning locations 

and ablation depths utilizing the simulated annealing approach may be analogous to those 
for the method that employs the least square fitting approach using analytical lens 
equations. 
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Stepl. Provide target profile matrix. The target profile matrix contains 
data that may represent any arbitrary target profile, or may be obtained by evaluating the 
lens equation at discrete points. 

Step 2. Define the basis function in the form of treatment matrix. The 
5 treatment matrix contains data that represent ablation profiles (i.e., sizes and shapes). The 
ablation profiles can be arbitrary as long as'they can be measured. 

Step 3. Define the positions for evaluations. The positions for 
evaluations may be selected, for instance, based on the target profile, the ablation profiles, 
or the like. 

10 Step 4. Define the merit function. Any desirable merit function may be 

specified. An example is a weighted least square fit with a standard deviation p\ specified 
for each data point. 

Step 5. Apply simulated annealing. The simulated annealing process as 
illustrated in Fig. 7 and described above is applied to fit the treatment matrix to the target 

15 profile matrix to solve for the treatment plan or table containing the amplitudes for a 
treatment partem which is in general a two-dimensional pattern. 

Step 6. Evaluate closeness of fit: The closeness of the fit is monitored by 
evaluating the merit function which can be a least square fit error function or any 
arbitrary function specified by the user. The evaluation occurs in step 26 in the flow chart 

20 shown in Fig. 7. 

Step 7. Iterate for best solution. As illustrated in Fig. 7, the iteration 
scheme involves random variation of the PIV and evaluation of the merit function (steps 
24-30), and adjusting the "temperature', i ; e., the size of incremental change to an element 
of the PIV (steps 32 and 34). 

25 Step 8. Determine actual pulse locations. The treatment plan or table 

obtained by solving the target profile fitting equation includes amplitudes of the pulses 
which, in conjunction with the treatment matrix, provides instructions to the laser to 
direct pulses with the sizes and depths at the desired locations to achieve the best fit to the 
target profile. 

30 Step 9. Randomize table. The pulses are randomized for treatment to 

reduce the periodicity caused by the finiteness i of the pulse number. 
II. Exemplary Scanning System 
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Fig. 8 A shows a block diagram of a scanning system 100 for implementing 
the method of the present invention to perform corneal sculpting. The scanning system 
100 includes a laser 102 for generating a laser beam and laser delivery system optics 104 
for delivering the laser beam to a target which is the cornea 106 of an eye. The laser 
delivery system optics 104 include various optical elements used to focus, modify, and 
direct the beam in a scanning mode to the cornea 106. The energy level of the laser beam 
is typically about 50-250 mJ/pulse. Commercially available systems having suitable laser 
and optics for scanning include, for example, the VISX TWENTY/TWENTY EXCIMER 
LASER S YSTEM, from VISX, Incorporated of Santa Clara, California; Lasersight, 
Orlando, Florida; Chiron Technolas, Germany; and Autonomous Technology, Orlando, 
Florida. 

Optionally, the system 100 may include a calibration tool 108 for 
determining the characteristics of the laser beam which can be used as feedback for 
calibrating the scanning system 100. In one embodiment, the calibration tool 108 
includes a reference-edge and a photodetector. The calibration tool 1 08 is placed at a 
. target location for sensing a laser beam directed onto the tool 108. The laser beam can be 
repeatedly re-directed between the tool i08 and the patient's cornea 1 06. As such, after 
determining the size, shape, and/or position of the beam using the tool 108, the laser beam 
can be applied at a known location on the cornea. A repetitive measurement of intensity 
and shape characteristics of the laser beam can be made, and a repetitive recalibration of 
the targeting of the laser beam can be achieved to ensure precise positional accuracy 
when ablating the cornea 106. Using the sensed information from the calibration tool 
108, an algorithm for calculating the locations and number of scanning spots (such as that 
according to the present invention) can be revised, thereby increasing the accuracy of the 
sculpting procedure. This calibration information can be used to adjust the ablation 
algorithm immediately before and/or during each ablation procedure in real time. A more 
thorough discussion of the calibration procedure is found in commonly assigned U.S. 
Patent Application entitled "Method and Apparatus for Determining Characteristics of a 
Laser Beam Spot," Serial No. 09/395,809 (Attorney Docket No. 18158-124), filed 
September 14, 1999, which is incorporated herein by reference in its entirety. 

The scanning system 100 further includes a controller 110 for controlling 
operation of the laser 102 and the delivery system optics 104, and all activities of the 
system 100. In one embodiment, the controller 110 includes a hard disk drive (memory 
112), a floppy disk drive, and a processor 114. The controller 110 executes system 
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control software, which is a computer program stored in a computer-readable medium 
such as the memory 1 12. The memory 1 12 is typically a hard disk drive, but may also be 
other kinds of memory. The computer program includes sets of instructions that dictate, 
for instance, spot ablation.pattern, spot ablation depth, depth per pulse, the energy level, 
5 spot size and shape, pulse rate, pulse location, and scanning pattern and sequence. It is 
understood that other computer programs stored on other memory devices including, for 
example, a floppy disk or another appropriate drive, may also be used to operate the 
controller 110. ; 

An interface 120 is provided between a user and the controller 110 

10 typically in the form of a display monitor for displaying information, and an input device 
such as a keyboard, a mouse, and/or a light pen to allow the user to communicate with the 
controller 1 10. ■' • 

The scanning process can be implemented using a computer program 
product that is executed by the controller 110. The computer program code may be 

15 written in any conventional computer readable programming language. Suitable program 
code is entered into a single file, or multiple files, using a conventional text editor, and 
stored or embodied in a computer usable medium, such as a memory system of the '• 
computer. If the entered code text is a high level language, the code isicpmpiled, and the 
resultant compiler code is then linked with an object code of precompiled library routines. 

20 To execute the linked, compiled object code, the system user invokes the object code, 
causing the computer system to load the code in memory. The processor 114 then reads 
arid executes the code to perform the tasks identified in the program. 

Fig. 8B shows an illustrative block diagram of the control structure of the 
system control software, computer program 130, according to a specific embodiment. 

25 Through the input device of the interface 120, a user enters a set of process parameters 
into a process control subroutine 1 32 in response to menus or screens displayed on the 
monitor. The process parameters include input data needed to operate the system 100 
including, for example, the laser beam energy and spot size, and all the input necessary to 
determine the pulse locations and numbers according to the algorithm described above. 

30 A scanning location determination subroutine 1 34 includes program code 

for accepting the process parameters from the process control subroutine 132, and for 
determining the scanning pulse locations and numbers according to the algorithm 
described aboye. The scanning location determination subroutine 134 produces laser 
control parameters and optics control. parameters used in operating the laser 102 and 
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delivery system optics 104. A laser control subroutine 136 includes program code for 
accepting laser control parameters to control operation of the laser 1 02. An optics control 
subroutine 1 38 includes program code for accepting optics control parameters to control 
operation of the system optics 104. A calibration control subroutine 140 includes 

5 program code for controlling the calibration tool 1 08 to determine the characteristics of 
the laser beam generated by the system 100 and to provide feedback data to calibrate the 
laser 102 and delivery system optics 104 for improved scanning accuracy and precision. 
The calibration data can also be used to adjust the ablation algorithm in the scanning 
location determination subroutine 134 to revise the calculation of scanning spot locations 

10 and numbers. 

The computer program 1 30 can be used to determine the pulse locations 
and numbers, and to calibrate the laser system 100 in real time. In this way, the target 
profile for a patient can be entered just prior to performing the corneal sculpting 
procedure. 
15 m. Examples 

The following experimental examples are used to illustrate the methods of 

the present invention. 

A. Least Square Fitting Approach Using Ana lytical Lens Equation 

The examples were undertaken using a scanning laser system having a 

20 Questek impulse excimer with a pre-objective galvo scanner. The pulse rate was 70 Hz. 
The treatment time is dependent on the amount of tissue removed per pulse (or tissue 
removed per second). For instance, a 6 mm spot will remove 9 times more material per 
pulse than a 2 mm spot. To keep the ablation rate constant, the 2 mm spot will need to 
increase the repetition rate by 9 as compared to the 6 mm spot. Decreasing the spot 

25 further will increase the number of pulses dramatically, as illustrated in the following 
table. 



Spot Size 


Volume ' , ' - 
rmm 3 l 


0.5 mm 


4.7x1 0" y 


1.0 mm 


1.9xl0- 8 


2.0 mm 


7.5x10 s 


6.0 mm 


6.8. 10"' 
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1. Myopia 

Fig. 9 is a graph comparing a fit using the method of the present invention 
to the theoretical lens equation. The number of basis, M, used in the fit is ; 27. Thelitis 
given a standard deviation, ct=0.25, for the first 1 50 points and a value four times larger 
5 for the last 1 50 points, cf=1 .0. Having ersmaljer for the inner portion of the curve gives 
the central 3 mm of a 6 mm ablation more importance than the outer portion. The 
^ difference in the quality of fit between the inner and the outer portions is evident, where 
the fit is much better in the inner portion. The outer portion shows a slight departure. 

One gross feature observed in this fit is the "bump" at the one-millimeter 
10 distance. This occurs because of the size of the spot (2 mm) is an integral multiple of the 
ablation zone (6 mm). One way to fix the problem is to increase the ablation zone 
.slightly such that the ablation zone is not an integer multiple of the spot. This will also 
taper the outer perimeter giving a smooth transition from the ablated area to the non- 
ablated area. Further, there are numerous small bumps near the center ablation zone in 
15 Fig; 9. These manifest themselves as concentric rings when the ablations are tested oh 
plastics. Increasing the number of basis functions will eliminate most of the rings. 

Shown in Fig. 10 is the same fit but with an increase from 27 to 47 basis 
functions. The ablation zone is not an integer multiple of the spot. The fit has the 
extended smooth transition between the ablated and non-ablafed zones. Most of the 
20 • undulations have disappeared. It should be noted, however, that care needs to be taken 
when increasing the number of basis functions, because too many basis function may 
cause the solution to require negative values making such solutions not physically 
obtainable. . 

Fig. 1 1 is an interferometer image of a -4 diopter ablation on a plastic 
25 using the simulation output from Fig. .10 (extended fit with 47 basis functions) and a spot 
size of approximately 2.1 mm. The salient 'feature is the circular over-ablation at 1 mm 
radius. ■ ■ - . 

Fig. 12 shows the characteristics of the simulated ablation profiles of a -4 
diopter ablation with 47 basis functions with variations in spot sizes from a 2 mm 
30 reference spot size. The curves above the 2 mm reference curve have 5, 10, 15 percent 
larger spot sizes (i.e., 2.1, mm, 2.2 mm, and 2.3 ntm), while the curves below have 5, 10, 
15 percent smaller spot sizes (i.e., 1.9 mm, 1.8 mm, and 1.7 mm). The dip or peak at 1 
mm becomes larger as the spot size deviates from nominal, and the under- or over- 
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ablation is non-linear. As a result, a change in spot size may remove less or more 
material and change the overall profile. 

2. Hyperopia 

Fig. 13 is an interferometer image of a hyperopic ablation profile of a +4 
diopter ablation with 33 basis functions using a 2 nun spot size for a 6 mm X 9.6 mm 
ablation zone. The method to generate hyperopic lenses is similar to that for myopia. 
The fitting requirement for hyperopia is less restrictive than for myopia. In myopia, both 
the ablation's central and perimeter contours (and everything in between) desirably 
follow the lens equation. In hyperopia, the ablation is a negative lens, so that the central 
portion has very few pulses and the perimeter has the majority of the pulses. The fitting 
in the central portion is much less stringent for hyperopia than for myopia. The pulses 
desirably are able to extend beyond the 6 mm treatment zone length in hyperopia. The 
outermost pulses have their centers at the 3 mm radius. For a 2 mm spot, it is desirable 
for the pulse to extend beyond the 3 mm radius going to about 3.5 mm to achieve a good 
fit. This extension produces an overall ablation diameter of about 9.6 mm. To reduce the 
ablation size, a smaller spot may be used- 
Figs. 14A and 14B are graphs of two hyperopic lenses employing the same 
basis functions used in the myopia examples and two different spot sizes. Fig. 14A uses a 
2 mm spot and Fig. 14B uses a 1.5 mm spot. The effective dimensions are 6x10 and 6x8 
mm ablations, respectively. A transition zone is observed from the 3 mm radius onwards. 
The spot size has a large influence on this transition zone. The larger the spot size, the 
larger the transition zone becomes. For a 6 mm optical zone, a 2 mm spot produces a 
transition zone of 4.8 mm radius, while a 1.5 mm spot produces a transition zone of 3.9 
mm radius. 

Fits to a hyperopic lens produce different characteristics from those 

obtained from the myopic fits. The hyperopic fits are relatively simple because the 

perimeters are allowed to "float" thereby producing a transition zone. As with myopia, an 

extended zone is desirable for the hyperopic fit. This extended zone is about 0.8 times the 

radius. This places the spots at 3.8 mm for a 2 mm spot. Adding the radius of the 2 mm 

spot gives the ablation zone radius of 4.8 mm. 

treatment zone . „ , . 0 >. 

= optical zone + 1.8 (spot size 12) 

radius 

Based on this equation, the treatment zone, diameter for a 1.5 mm spot is 8.7 mm. 
3. Hyperopic Cylinder 
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Using the appropriate lens equation and basis functions, a hyperopic 
cylinder ablation profile is generated for a +4 diopter ablation with 33 basis functions 
using a 2 mm spot size for a 6 mm X 9.6 mm ablation zone, as shown in Fig, 15. 

Employing the method of the invention, the fit returns the amplitudes of 

5 the basis functions. Stacking the pulses vertically will give the correct depth but it will 

not produce the lateral dimensions needed. As described earlier, offsetting the pulses will 

produce the trench. The boundary conditions are given by the rectangular perimeter of 

the hyperopic cylinder. For example, a 5x9 mm hyperopic cylinder has a 9 mm width 

along the lens direction (similar to plain hyperopia) and has a 5 mm cut along the 

1 0 perpendicular directions. The 5 mm is the length of the trench not including the transition 

zone. The transition zone is the overlap-offset used to generate the basis "trenches." The 

slop of the transition zone is given by 

s j 0 p e _ depth per pulse 
• * offset 

The offset is constrained by the total depth, the spot size, and trench 

15 length. The relationship among the total depth, spot size, and offset is given by : 

. spot size 

step size = — r 

depth/depth per pulse 

The transition zone is one spot size in width. The number of the total 

pulses for this basis trench is not simply the depth/depth per pulse, since it also depends 

on the trench length. Instead, the number of pulses is 

_ n _ . length of trench depth 

20 no. of pulses = — — — -*— — — 

spotsize depth per pulse 

For myopia and hyperopia, the periodicity is reduced considerably by 
adding a random starting phase to each ring. This approach is not available to ablation 
profiles with cylindrical symmetry such as hyperopic and myopic cylinders. Starting the 
trench with a repetitious position will produce a deeper pattern at that position with a 

25 feathering effect at the ending positions. The feathering effect occurs because the number 
of pulses per basis is different making the offsets different. A very small random phase 
or position is desirably added to reduce the periodicity of the scans. For instance, a small 
random edges position is added to the locations to blend the spot. 

Generating a myopic cylinder ablation profile is similar to generating a 

30 hyperopic cylinder profile. The lens equations are different but the same basis functions 
are used. Fitting the basis functions to the lens equation for a myopic cylinder is more 
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stringent than for a hyperopic cylinder, however, because the fit is preferably good from 
the center to the perimeter with a small allowance for the transition zone. The fit in the 
center region for a hyperopic cylinder is less rigorous. 

In addition, care needs to be taken when the treatment ring approaches the 
origin. Analytically, the basis functions for the myopic cylinder are equivalent to the 
ones used for hyperopic cylinder. The difference is the hyperopic cylinder treatment does 
not approach the origin. At the origin, the basis functions overlap onthemselves. The 
overlap is 

y. O) = J(spotsize/2) 2 -(*-x o .) 2 

+ J(spot size / 2) 2 - (x + x oi ) 2 

where spot size /2 <(x oi + x). This "overlap" is part of the algorithm and is used if the fit 
goes from 0 to the treatment zone radius, to account for the situation where the spot 
passes through the original ([0]). 

B. Simulated Annealing Approach 

An example employs the simulated annealing approach for fitting the 
ablation of a -2 diopter elliptical profile which has a major axis length of 16 mm and a 
minor axis length of 12 mm, and a maximum depth of about 3 1 um. Due to the bilateral 
symmetry, only one quarter of the surface is used in the fitting calculations. Any overlaps 
caused by large pulses near the axes of symmetry are reflected back relative to the axes of 
symmetry to ensure that the calculations are correct. The pulse shapes are circular having 
diameters of 5.5 mm, 5 mm, 4 mm, 3 mm, 2 mm, and 1 mm at evenly spaced (100 um) 
intervals. There are about 2000 total pulse basis functions. The merit function in this 
case is the root mean square error between the target shape and the current solution. 

Figs. 16a- 16c show the simulated annealing process at different stages of 
the iteration to achieve the best fit. Each figure shows (1) a working solution obtained 
from the current best guess of the treatment plan or table, (2) the target profile, and (3) 
delta representing the difference between the working solution and the target profile. At 
the initial guess with zero iteration shown in Fig. 16a, the errors are quite large. The 
errors are reduced to about 2 um or less after 1 0,000 iterations as shown in Fig. 1 6b. At 
1,000,000 iterations in Fig. 16c, the errors are less than 1 um and well distributed over the 
surface. The contour plots of the target profile and the calculated profile are illustrated in 
Figs. 17a and 17b, respectively. About 50% of the calculated profile falls within 0.1 urn 
of the target profile, and about 80% of the calculated profile falls within 0.25 um of the 
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T 

target profile. The solution is used by a Tencor scanner to produce the ablation profile 
shown in Fig. 18. 

It is to be understood that the above description is intended to be 
illustrative and not restrictive. Many embodiments will be apparent to those of skill in the 
5 art upon reviewing the above description. For instance, some of the above equations are 
developed using a fixed spot shape and size. In alternate embodiments, the spot shape 
and/or spot size may be varied. Variation of spot size and shape is described, for 
instance, in commonly assigned European Patent No v 0628298 for "METHOD AND 
SYSTEM FOR LASER TREATMENT OF REFRACTIVE ERRORS USING OFFSET 

10 IMAGING", filed May 5, 1994 (Application No. 94303256.5) and issued April 1, 1998, 
which is incorporated herein by reference in its entirety. Different merit functions can be 
prescribed. Other numerical approaches may be used to solve the target profile fitting 
equation. The desired lens profile may be an elliptical profile, a hyperopie elliptical 
profile, a myopic elliptical profile, a circular profile, a linear profile, an asymmetric 

15 profile, or an arbitrary two-dimensional lens profile. The scope of the invention should, 
therefore, be determined not with reference to the above description, but instead should 
be determined with reference to the appended claims along with their full scope of 
equivalents. 
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WHAT IS CLAIMED IS: 

1 1 . A method of generating a treatment table for ablating tissue using a 

2 scanning laser beam for generating scanning spots over a treatment region larger in area 

3 than the scanning spots, the method comprising: 

4 providing a target function representing a desired lens profile for ablating 

5 the tissue by scanning spots of the laser beam on the tissue; 

6 providing a basis function representing a treatment profile produced by 

7 scanning with overlapping scanning spots of the laser beam in a treatment pattern; and 

3 fitting the target function with the basis function to obtain a treatment table 

9 including scanning spot locations and characteristics for the overlapping scanning spots 

10 of the laser beam. 

1 2. The method of claim 1 wherein the basis function is a two- 

2 dimensional function representing a two-dimensional section of a three-dimensional 

3 treatment profile which has symmetry with respect to the two-dimensional section 

4 extending along the treatment pattern. 

1 3, The method of claim 2 wherein the treatment pattern is generally 

2 linear or generally circular. 

1 , 4. The method of claim 1 wherein the target function is a two- 

2 dimensional function representing a two-dimensional section Of a three-dimensional lens 

3 profile which has symmetry with respect to tne two-dimensional section extending along 

4 the treatment pattern. 

1 5. The method of claim 4 wherein the target function represents an 

2 ablation depth as a function of a distance from an optical axis of a cornea. 

1 6. The method of claim 1 wherein fitting the target function with the , 

2 basis function includes fitting at N discrete evaluation points. 

1 7. The method of claim 6 wherein the basis functibn includes M 

2 discrete basis functions representing M overlapping scanning spots. 

j 8. The method of claim 7 wherein the M discrete basis functions 

2 represent M overlapping scanning spots across a treatment zone length representing the 

3 length across a generally two-dimensional section which is oriented normal across a 
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4 generally straight treatment pattern or which is oriented radially across a generally 

5 circular treatment pattern. 

1 9. The method of claim 8 wherein the scanning spots are generally 

2 circular and have a generally uniform energy profile. 

1 10. The method of claim 9 wherein 

2 (A) for a treatment profile having a generally uniform two-dimensional 

3 section oriented normal across a generally straight treatment pattern, the discrete basis 

4 functions represent the two-dimensional section as 



5 Xfa) =>,(*,) = Jis/2) 2 -(xj-x^Y or 

6 (B) for a treatment profile having a generally uniform two-dimensional 

7 section oriented radially across a generally circular treatment pattern, the discrete basis 

8 functions represent the two-dimensional section as 

V + *o,-(s/2) :N 



^,(*,) = 3 (*,) = cos" 1 



2x o; Xj 



10 where 

11 s is the diameter of the scanning spot; 

12 . j= 1 N; 

13 xj is a reference ^-coordinate for the two-dimensional section measured 

14 from an optical axis of the cornea of a j tt? evaluation point for the center of the scanning 

15 spot; 

16 xoi is an ^-coordinate for a center of an i ,h scanning spot; 

17 (xoi - s/2) < Xj < (xoi + s/2); 

18 yi(xj) is a depth of the i th basis function for the generally straight treatment 

19 pattern; and 

20 9i(Xj) is a coverage angle of the i ,h basis function for the generally circular 

21 treatment pattern. 

1 11. The method of claim 10 wherein x 0l is specified for M number of 

2 equally spaced scanning spots as = i * [(L - s + e) / M], 

3 . where 

4 /. is the treatment zone length; 

5 e is an extended zone; and 

6 / = 1,...,M. ' 
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1 12. The method of claim 11 wherein e is set to about 0.1 to about 0.5 

2 mm. 

\ 13. The method of claim 7 wherein M is equal to about 7 to about 97. 

1 14. The method of claim 7 further comprising refitting the target 

2 function with the basis function by varying the number of scanning spots M to iterate for 

3 a best fit. 

1 15. The method of claim 6 wherein the target function is: 

2 (A) for myopia and myopic cylinder, 



5 



11 



n-X + R t D J J 
4 (B) for hyperopia and hyperopic cylinder, 



-x 2 ; + C or 



' — X ; or 



6 (C) for photo therapeutic keratectomy, 

.7 f(xj)=d; 

8 where 

9 0<Xj<_(L -shift); 

10 y = o,l,...,N-lv 



12 x,- is an ^-coordinate measured from an optical axis of the cornea of the j* 

13 evaluation point for the center of the scanning spot; 

14 s is the diameter of the scanning spot; 

15 R, is the anterior radius of curvature of the cornea in meters; 

16 R 2 is the final anterior radius of curvature of the comea in meters; 

17 n = 1.377 is the index of refraction of the comea; 

18 D is the lens power of the scanning spot in diopters; 

19 L is the treatment zone length representing the length across a generally 

20 uniform section which is oriented normal across a generally straight treatment pattern for 

21 myopic or hyperopic cylinders, or which is oriented radially across a generally circular 

22 treatment pattern for myopia or hyperopia; 
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23 ■ shift is the amount of emphasis shift; and 

24 d is a constant depth. 

1 16. The method of claim 15 wherein the shift is about 0 to about 0.2. 

1 17. The method of claim 15 wherein xj=j *[ (L - shift) I N]. 

1 18. The method of claim 15 wherein the basis function includes M 

2 discrete basis functions representing M overlapping scanning spots, and wherein fitting 

3 the target function with the basis function comprises solving the following equation for 

4 coefficients a, representing treatment depth for the i ,h scanning spot: 

• - <=i 

6 where 

7 Xi(xj) is the i th basis -function; and 

8 i= 1,...,M. 

1 The method of claim 6 wherein fitting the target function and the 

2 basis function comprises specifying a deviation for each of the N discrete evaluation 

3 points. 

t 20. The method of claim 19 further comprising refitting the target 

2 function with the basis function by varying the deviations to iterate for a best fit. 

1 21. The method of claim 1 wherein fitting the target function and the 

2 basis function comprises evaluating closeness of the fit and repeating the fitting step if the 

3 closeness does not fall within a target closeness. 

1 22. The method of claim 1 wherein the target function and the basis 

2 function are fitted using a least square fit. 

1 .23. The method of claim 1 further comprising randomizing the 

2 scanning spot locations of the treatment table to produce a random scanning order; 

1 24. The method of claim 1 further comprising refitting the target 

2 function with the basis function by varying the size of at least one of the scanning spots to 

3 iterate for abest fit. 
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1 25. The method of claim 1 wherein the scanning spot characteristics of 

2 a scanning spot at a scanning spot location include shape, size, and depth of the scanning 

3 spot at the scanning location. . 

1 26. The method of claim 1 wherein the scanning spots have different 

2 sizes. . 

1 27. The method of claim 1 further comprising specifying the treatment 

2 pattern for scanning with overlapping scanning spots of the laser beam. 

1 28. The method of claim 1 wherein the target function and the basis 

2 function are fitted using a simulated annealing process. 

1 29. The method of claim 1 further comprising specifying a merit 

2 function representing an error of fit between the target function and the basis function; 

3 and minimizing the merit function. 

1 '30. The method of claim 1 further comprising specifying a merit 

2 function representing an error of fit between the target function and the basis function; 

3 monitoring a total number of the scanning spots in the treatment table; and minimizing 

4 the merit function and the total number of the scanning spots in the treatment table. 

1 31, The method of claim 1 further comprising refitting the target 

2 function with the basis function by selecting a scanning spot location and varying the 

3 characteristics of the scanning spot at the selected location to iterate for a best fit. 

1 32. A method of generating a treatment table for ablating tissue using a 

2 scanning laser beam for generating scanning spots over a treatment region larger in area 

3 than the scanning spots, the method comprising: 

4 providing a lens function representing a desired lens profile for ablating 

5 the tissue by scanning spots of the laser beam on the tissue; 

6 providing a basis function representing a treatment profile produced by the 

7 overlapping scanning spots along a treatment path, the basis function representing a 

8 section oriented across, the treatment path; and 
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9 fitting the lens function with the basis function to obtain a treatment table 

10 including scanning spot locations and characteristics for the overlapping scanning spots 

11 of the laser beam. 

1 33. The method of claim 32 wherein the scanning spots are generally 

2 circular and have a generally uniform energy profile, and the basis function includes M 

3 discrete basis functions representing M overlapping scanning spots. 

1 34. The method of claim 33 wherein the treatment profile is 

2 symmetrical with respect to an axis of symmetry, and the discrete basis functions are 

V+4-(s/2) jX 



■ 9.{x) - cos" 



2-x oi x 



4 where . 

5 ... s is the diameter of the scanning spot; 

6 . x is an jc-coordinate measured from the axis of symmetry; 

7 xoi is an x-coordinate for a center of an i ,h scanning spot; , 

8 (x 0 i - s/2) <x< (x 0i + s/2); and 

9. 6i(x) is a coverage angle of the i lh basis function. 

1 35. The method of claim 34 wherein .xp, is specified for M number of 

2 equally spaced scanning spots as: 

3 x 0i = i*{(L-s + e)/M}, 

4 where 

5 L is the treatment zone length of the section oriented radially across the 

6 treatment profile; 

.7. ' e. is an extended zone; and 

.8 '. i= 1,...,M. / ■ ' , ;'. 

1 36. The method of claim 34 wherein fitting the leiis function with the 

2 basis function comprises solving the following equation for coefficients a,- representing 

3 treatment depth for the i ,h scanning spot: 



f(x) = f j a l X i (x) 



5 where . 

6 f(x) is the lens function; and 



38 



WO 01/67978 



PCT/US01/08337 



7 i= 1„..,M. 

1 37. The method of claim 36 wherein the lens function is: 

2 (A) for myopia, 



10. 



4 (B ) for hyperop i a, 



6 (C) for phototherapeutic keratectomy, 

7 

8 where 

9 0 <x <_(L - shift); 



-x 2 or 



C = jR?-s z /4 + 



11 5 is the diameter of the scanning spot; 

12 i? ; is the anterior radius of curvature of the cornea in meters; 

13 i? 2 is the final anterior radius of curvature of the cornea in meters; 

14 n = 1.377 is the index of refraction of the cornea; 

15 D is the lens power of the scanning spot in diopters; 

16 L is the treatment zone length; 

17 shift is the amount of emphasis shift; and 

18 d is a constant depth. 

1 38. The method of claim 36 further comprising dividing the depth (a,) 

2 for the i th scanning spot by a depth per pulse of the laser beam to obtain a number of 

3 pulses per an i ft treatment ring for the i th scanning spot; and dividing the number of pulses 

4 per treatment ring by 2n to obtain an angular spacing between pulses for the i m treatment 

5 ring. 

\ 39. The method of claim 32 wherein the scanning spots have a fixed 

2 spot size and a fixed spot shape. 
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1 40. The method of claim 32 wherein at least one of the spot size and 

2 spot shape of the scanning spot is variable. 

1 41. A method for fitting a three-dimensional target profile, the method 

2 comprising: 

3 providing a two-dimensional basis function including overlapping portions 

4 to represent a three-dimensional profile which has symmetry with respect to a two- 

5 dimensional section extending along a treatment pattern; and 

6 fitting the three-dimensional target profile with the two-dimensional basis 
.7 function to obtain a distribution of the overlapping portions. 

1 42. The method of claim 41 wherein the three-dimensional profile has 

2 symmetry with respect to a two-dimensional section oriented radially from an axis of 

3 symmetry and extending in a generally circular treatment pattern around the axis. 

.1 . 43. The method of claim 42 wherein the overlapping portions are 

2 generally circular, arid the two-dimensional basis function comprises discrete basis 

3 functions each representing a coverage angle of one of the overlapping portions as a 

4 function of a distance from the axis of symmetry. 

1 44. The method of claim 41 wherein the three-dimensional profile has 

2 symmetry with respect to a two-dimensional section oriented normal across a generally 

3 straight treatment pattern. 

1 45. The method of claim 44 wherein the overlapping portions are 

2 generally circular, and the two-dimensional basis function comprises discrete basis 

3 functions each representing a depth of one of the overlapping portions as a function of a 

4 distance from the axis of symmetry; 

1 46. A system for ablating tissue, the system comprising: 

2 a laser for generating a laser beam; , 

3 a delivery device for delivering the laser beam to a tissue; 

4 a controller configured to control the laser and the delivery device; and : 

5 a memory, coupled to the controller, comprising a computer-readable 

6 medium having a computer-readable program embodied therein for directing operation of 
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7 the system, the computer-readable program including a first set of instructions for 

8 generating a treatment table including scanning spot' locations and characteristics for 

9 ablating the tissue over a treatment region larger in area than the spot size of the laser 

1 0 beam to achieve a desired lens profile for ablating the tissue, a second set of instructions 

1 1 for controlling the laser to generate the laser beam, and a third set of instructions for 

12 controlling the delivery device to deliver the laser beam to the tissue according to the 

13 treatment table. 

1 47. The system of claim 46 wherein the first set of instructions of the 

2 computer-readable program includes: 

3 a first subset of instructions for providing a target function representing the 

4 desired lens profile for ablating the tissue by scanning spots of the laser beam on the 

5 tissue; 

6 a second subset of instructions for providing a basis function representing 

7 a treatment profile produced by the overlapping scanning spots in a treatment pattern; and 

8 a third subset of instructions for fitting the target function with the basis 

9 function to obtain the treatment table including the scanning spot .locations and 

1 0 characteristics for the overlapping scanning spots of the laser beam. 

1 ■ 48. The system of claim 47 wherein the second subset of instructions 

2 provide a basis function which is a two-dimensional function representing a two- 

3 dimensional section of a three-dimensional treatment profile having symmetry 

4 respect to the two-dimensional section extending along the treatment pattern. 

1 49. The system of claim 47 wherein the first set of instructions of the 

2 computer-readable program includes a fourth subset of instructions for refitting the target 

3 function with the basis function by varying the spot size of the laser beam to iterate for a 

4 best fit. 

1 50. The system of claim 47 wherein the first set of instructions of the 

2 computer-readable program includes a fifth subset of instructions for evaluating closeness 

3 of the fit and repeating the fitting step if the closeness does not fall within a target 

4 closeness. 
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1 . 51. The system of claim 47 wherein the first set of instructions of the 

2 computer-readable program includes a sixth subset of instructions for randomizing the 

3 scanning spot locations for the treatment table to produce a random scanning order. 

1 52. The system of claim 47 wherein the first set of instructions of the 

2 computer-readable program includes a seventh subset of specifying the treatment pattern 

3 for scanning with overlapping scanning spots of the laser beam; 

1 53. The system of claim 47 wherein the scanning spot characteristics of * 

2 a scanning spot at a scanning location include shape, size, and depth of the scanning spot 

3 at the scanning location. 

1 54. The system of claim 47 wherein the desired lens profile is selected 

2 from the group consisting of an elliptical profile, a hyperopic elliptical profile, a myopic 

3 elliptical profile, a circular profile, and a linear profile. 

1 55. The system of claim 47 wherein the desired lens profile. is 

2 asymmetric. 

1 56. The system of claim 47 wherein the desired lens profile comprises 

2 an arbitrary two-dimensional lens profile. 
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